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ABSTRACT 

The observed spectra of blazars, their intrinsic emission, and the underlying populations of radiat- 
ing particles are intimately related. The use of these sources as probes of the extragalactic infrared 
background, a prospect propelled by recent advances in TeV-band telescopes, soon to be augmented 
by observations by NASA’s upcoming Gamma-Ray Large Area Space Telescope (GLAST), has been a 
topic of great recent interest. Here, it is demonstrated that if particles in blazar jets are accelerated at 
relativistic shocks, then 7 -ray spectra with indices less than 1.5 can be produced. This, in turn, loosens 
the upper limits on the near infrared extragalactic background radiation previously proposed. We also 
show evidence hinting that TeV blazars with flatter spectra have higher intrinsic TeV 7 -ray luminosities 
and we indicate that there may be a correlation of flatness and luminosity with redshift. 

Subject headings: gamma-rays: theory, (galaxies:) BL Lacertae objects: general, (cosmology:) diffuse 
radiation 

1. INTRODUCTION 

The comparison of theoretical models for the 7 -ray spectra 
of blazars with observations is the standard approach used to 
understand the physical processes leading to their high energy 
emission. In the case of TeV components, the radiation models 
generally involve the synchrotron self-Compton (SSC) process. 

Determining the intrinsic emission spectra of such TeV sources 
requires that one account for the energy and redshift dependent 
absorption of 7 -rays from these sources through 77 — > e+e“ in- 
teractions with intergalactic photon backgrounds produced by 
stellar and dust emission. The cumulative background radiation 
seen at redshift z = 0 is commonly referred to as the extragalac- 
tic background light (EBL). Various calculations of extragalac- 
tic 7 -ray absorption have been given in the recent literature and 
they are discussed along with the latest calculations in the paper 
of Stecker, Malkan & Scully (2006, hereafter SMS06). 

The results given by SMS06 are based on two galaxy evo- 
lution models, viz. a baseline model (B) and a fast evolution 
model (FE). The spectral energy distributions of the extragalac- 
tic background light for these models are shown in Figure 1 . The 
FE model is favored by recent Spitzer observations (Le Floc’h 
et al. 2005, Perez-Gonzalez et al. 2005). It provides a bet- 
ter description of the deep Spitzer number counts at 70 and 
160/im than the B model. However, GALEX observations indi- 
cate that the evolution of UV radiation for 0 < z < 1 may be 
somewhat slower and more consistent with the B model within 
errors (Schiminovich et al. 2005). And the 24/xm Spitzer source 
counts are closer to the B model than the FE model. 

Probing the EBL using blazar observations is contingent upon 
an accurate understanding of their emission spectra. Aharonian 
et al. (2006) have argued that intrinsic blazar spectra must have 
spectral indices T s > 1.5. They use this assumption, together 
with HESS observations of the source 1 ES 1101-232, to place 
an upper limit on the EBL of 14 nWm _ 2 sr _1 at a near infrared 


wavelength of 1.5 /im, corresponding to a frequency of 2 x 10 14 
Hz. As can be seen from both Figure 1 and Table 1 , this value 
is consistent with model B, but not with the model FE, which 
is favored by the Spitzer observations. It is therefore important 
for exploring both galaxy evolution and blazar physics that we 
reexamine the assumption T s > 1.5 made by Aharonian et al. 
(2006). This assumption has been questioned by Katarz’nski et 
al. (2006) in a different context, but we will examine it here in 
the light of the physics of shock acceleration, which provides in- 
sights into the distributions of underlying particle populations. 



Fig. 1 . — Spectral energy distribution of the EBL, taken from SMS06. 
Data points with error bars depict measurements, triangles show lower lim- 
its from source number counts, and the inverted triangle shows an upper 
limit from Stecker and De Jager (1997). The upper and lower solid curves 
depict FE (fast evolution) and B (baseline) model predictions. The dotted 
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Table 1 

Blazar Spectral Indices in the 0.2 — 2 TeV Energy Range and Isotropic Luminosities at 1 TeV 


Source 

z 

T obs 

r s (FE -» B) 

T e (FE -» B) 

£(1 TeV) [10 36 WJ 

1ES 2344+514 

0.044 

3.0 

2.5 -» 2.6 

4.0 -► 4.2 

2.9 

Mrk 180 

0.045 

3.3 

2.9 -> 3.0 

4.8 -> 5.0 

1.2 

1ES1959+650 

0.047 

2.7 

2.3 -> 2.4 

3.6 -> 3.8 

5.4 

PKS 2005-489 

0.071 

4.0 

3.4 -> 3.5 

5.8 -> 6.0 

8.6 

PKS 2155-304 

0.117 

3.3 

2.2 -* 2.4 

3.4 -> 3.8 

420 

H 2356-309 

0.165 

3.1 

1.5 -> 1.9 

2.0 -> 2.8 

200 

1ES 1218+30 

0.182 

3.0 

1.2 -> 1.6 

1.4 -> 2.2 

310 

1ES 1101-232 

0.186 

2.9 

1.0 -» 1.5 

1.0 -> 2.0 

230 

1ES 0347-121 

0.188 

3.1 

1.2 -> 1.7 

1.4 -> 2.4 

1200 

1ES 1101+496 

0.212 

4.0 

1.8 -» 2.4 

2.6 -> 3.8 

930 


lines show the extensions of the models into the optical-UV (from SMS06). 

2. intergalactic ABSORPTION 

The intergalactic 7-ray absorption coefficient (i.e. optical 
depth), t(E, z ), increases monotonically with energy and there- 
fore leads to a steepening of the intrinsic source spectra as 
observed at Earth. SMS06 give a useful parametric form for 
t(E,z) with the corrected parameters given in the erratum 
(Stecker, Malkan & Scully 2007). For sources at redshifts be- 
tween 0.05 and 0.4, Stecker &: Scully (2006, hereafter SS06) have 
shown that this steepening results in a well-defined increase in 
the spectral index of a source with an approximate power-law 
spectrum in the 0.2 - 2 TeV energy range with spectral index 
T obs . This increase is a linear function in redshift z of the form 
A T = C 4- Dz , where the parameters C and D are constants. 
The overall normalization of the source spectrum is also reduced 
by an amount equal to exp{— (A + Bz )}, again where A and B 
are constants. The values of A, B, C, and D are given for the B 
and FE models in SS06. SS06 have used this relation to calcu- 
late the intrinsic 0.2 - 2 TeV power-law 7-ray spectra of sources 
having known redshifts in the 0.05 - 0.4 redshift range for both 
the B and FE models of EBL evolution. A version of Table 2 
of SS06 giving values for intrinsic spectral index of the source 
T s is shown in Table 1. Table 1 also shows the respective in- 
dices T e = 2 T s — 1 of the electron distributions in the sources 
under the assumption that the 7-rays are produced by inverse 
Compton interactions in the Thomson regime. 

Using the formula derived in SS06, we can estimate the in- 
trinsic “isotropic luminosity.” 1 The isotropic luminosity of the 
blazar sources listed in Table 1 is obtained from the formula 

C ~ 47T^|(1 + z) r °- 2 Fo[d( Z )] 2 e ( A+Bz > (1) 

where d is the luminosity distance to the source, and F 0 is its 
observed differential energy flux at 1 TeV. The other factors 
in the equation give the k-correction for the deabsorbed source 
spectrum and the normalization correction factor for absorp- 
tion given in SS06. The observational references for the sources 
listed in Table 1 can be found in SS06 except for the new obser- 
vations of 1ES 2344+514 (Albert et al. 2007a), 1ES 1959 + 650 
(Albert et al. 2006), 1ES 0347-121 (Aharonian et al. 2007) and 


1ES 1101+496 (Albert et al. 2007b). The source PG 1553+113 
at z = 0.36 (included in Table 2 in SS06) is not listed here 
because the observations are in the energy range 0.09-0.6 TeV 
and are therefore below the operative energy range for applying 
the analytic approximation given in SS06. The blazars Mrk 421 
and Mrk 501 are not included because their redshifts are signif- 
icantly less than 0.05. However, these blazars are analysed by 
Konopelko et al. (2003). 

The numbers given in the last column of Table 1 are de- 
rived for the fast evolution (FE) model. One may note that 
there appears to be a trend toward blazars having flatter in- 
trinsic TeV spectra and higher isotropic luminosities at higher 
redshifts. However, one must be careful of selection effects. The 
TeV photon fluxes of these sources as observed by HESS and 
MAGIC only cover a dynamic range of a factor of ~20. There- 
fore, only brighter sources can be observed at higher redshifts. 
This is because of both diminution of flux with distance, and in- 
tergalactic absorption (Stecker, de Jager & Salamon 1992). The 
observed C{z) trend is naturally expected in a limited popula- 
tion sample spanning a range of redshifts if the TeV-band fluxes 
are pegged near an instrumental sensitivity threshold. A more 
powerful handle on the intrinsic spectra and luminosities of these 
sources will be afforded by the upcoming GLAST 7-ray mission, 
with its capability for detecting many blazars at energies below 
200 GeV. 

3. PARTICLE DISTRIBUTIONS FROM SHOCK ACCELERATION 

As discussed above, inferences of source spectra from specific 
blazars are contingent upon the particular choice of an EBL 
model. Hence there is significant uncertainty in deductions of 
the underlying distribution of emitting electrons in the case of 
SSC models, or protons in the case of hadronic models. It is 
a goal of this presentation to provide a cohesive connection be- 
tween the particle distributions in blazars, the resulting emission 
spectra in the TeV band, and the spectrum of the EBL. 

The rapid variability seen in TeV flares drives the prevailing 
picture for the blazar source environment, one of a compact, 
relativistic jet that is structured on small spatial scales that 
are unresolvable by present 7-ray telescopes. Turbulence in the 
supersonic outflow in these jets naturally generates relativis- 
tic shocks, and these form the principal sites for acceleration of 
electrons and ions to the ultrarelativistic energies implied by the 


1 We define isotropic here as if the source had an apparent isotropic luminosity even though blazars are highly beamed and their flux (and hence their 

apparent luminosity) is dramatically enhanced by relativistic Doppler boosting. This is similar to the nomenclature used for 7 -ray bursts. The quantity 

C is equal to 4ni/F u given at hi/ = 1 TeV in units of 10 36 W (lCr 3 erg s -1 ). 
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TeV 7-ray observations. Within the context of this relativistic, 
diffusive shock acceleration mechanism, numerical simulations 
are used here to derive expectations for the energy distributions 
of particles accelerated in blazar jets. 

Diffusive acceleration at relativistic shocks is less well studied 
than that for nonrelativistic flows, yet it is the most applicable 
process for extreme objects such as pulsar winds, jets in ac- 
tive galactic nuclei, and 7-ray bursts. Early work on relativistic 
shocks was mostly analytic in the test-particle approximation 
( e.g ., Peacock 1981, Kirk & Schneider 1987, Heavens & Drury 
1988), where the accelerated particles do not contribute sig- 
nificantly to the global hydrodynamic structure of the shock. 
A key characteristic that distinguishes relativistic shocks from 
their non-relativistic counterparts is their inherent anisotropy 
due to rapid convection of particles through and away down- 
stream of the shock. This renders analytic approaches more dif- 
ficult for ultrarelativistic upstream flows, though advances can 
be made in special cases, such as the limit of extremely small 
angle scattering (pitch angle diffusion) (e.g. Kirk & Schneider 
1987; Kirk et al. 2000). Accordingly, complementary Monte 
Carlo techniques have been employed for relativistic shocks by 
a number of authors, including test-particle analyses for steady- 
state shocks of parallel and oblique magnetic fields by Ellison, 
Jones & Reynolds (1990), Ostrowski (1991), Bednarz h Os- 
trowski (1998), Baring (1999), and Ellison & Double (2004). It 
is such a simulational approach that is employed here to illus- 
trate key spectral characteristics for particles accelerated to high 
energies at relativistic shocks that are germane to the blazar 
emission-EBL attenuation problem. For a recent discussion of 
relativistic shock acceleration, see Baring (2004). 

The simulation used here to calculate diffusive acceleration in 
relativistic shocks is a Monte Carlo technique that has been 
employed extensively in supernova remnant and heliospheric 
contexts, and is described in detail in papers by Ellison, Jones 
& Reynolds (1990), Jones & Ellison (1991), Baring, Ellison h 
Jones (1994) and Ellison, Baring, Sz Jones (1996). It is concep- 
tually similar to Bell’s (1978) test particle approach to diffu- 
sive shock acceleration. Particles injected upstream gyrate in a 
laminar electromagnetic field, and particle trajectories are de- 
termined by solving a relativistic Lorentz force equation in the 
frame of the shock. Because the shock is moving with a velocity 
u relative to the plasma rest frame, there will, in general, be 
a u x B electric field in addition to the bulk magnetic field. 
Alfven wave turbulence is modeled by using a phenomenologi- 
cal description of ion scattering in the rest frame of the plasma. 
The scattering allows particles to diffuse spatially along mag- 
netic field lines, and to varying extent, across them as well. 
The scatterings are also assumed to be elastic, an assumption 
that is valid so long as the flow speed far exceeds the Alfven 
speed. Hence, contributions from stochastic second-order Fermi 
acceleration, where the scattering centers move with the Alfven 
waves, are generally neglected. The diffusion permits a minority 
of particles to cross the shock plane numerous times, gaining en- 
ergy with each crossing via the shock drift and first-order Fermi 
processes. 

A continuum of scattering angles, between large-angle or 
small- angle cases, can be modeled by the simulation. Denot- 
ing local fluid frame quantities by a subscript / , the time, Stf , 
between scatterings is determined by the mean free path, A / , 
the speed of the particle, Vf , and the maximum scattering an- 
gle, 0 S catt , as derived in Ellison, Jones, & Reynolds (1990); for 
small angles, it is given by Stf = A/0g Catt /(6t;/) , a formula that 
generally holds to within 10% for 0 sca tt < 80°. Here A / is 


proportional to a power of the particle momentum p (e.g., see 
Ellison et al., 1990; Giacalone et al., 1992 for microphysical jus- 
tifications for this choice), and for simplicity we assume that it 
scales as the particle gyroradius, r g , i.e. A/ = rjr g oc p. 

The parameter rj in the model is a measure of the level of tur- 
bulence present in the system, coupling directly to the amount of 
cross-field diffusion such that rj — 1 corresponds to the isotropic 
Bohm diffusion limit. It can be related to parallel ( /sy = X/v/3 ) 
and perpendicular ( k j_ ) spatial diffusion coefficients through 
the relation k_ l/k\\ = 1/(1+ rj 2 ) (see Forman, Jokipii & Owens, 
1974, Ellison et al., 1995, or Jokipii, 1987). In the parallel shocks 
considered here, where the B field is directed along the shock 
normal, rj has only limited impact on the resulting energy spec- 
trum, principally determining the diffusive scale normal to the 
shock. However, in oblique relativistic shocks, for which the 
field is inclined to the shock normal, the diffusive transport of 
particles across the field (and hence across the shock) becomes 
critical to retention of them in the acceleration process. Accord- 
ingly, for such systems, the interplay between the field angle and 
the value of rj controls the spectral index of the particle distri- 
bution (e.g. see Ellison and Double, 2004; Baring, 2004). 



Fig. 2. — Particle distribution functions dN/dp from parallel (©bi = 
0 ° ), relativistic shocks of upstr earn- to-downstr earn velocity compression 
ratio r = ui/u 2 = 3, as obtained from a Monte Carlo simulation of par- 
ticle diffusion and gyrational transport. Three shock speeds u\ = ftc 
are depicted, namely ft = 0.9487 , ft = 0.995 and ft = 0.99944 , corre- 
sponding to the labels 71 ft = 3 , 71 ft. = 10 and 71 ft = 30 , respectively, 
on the heavyweight (blue) histograms. Scattering off hydromagnetic tur- 
bulence is modeled by randomly deflecting particle momenta by an angle 
within a cone, of half-angle 0 S catt , whose axis coincides with the particle 
momentum prior to scattering. The heavyweight (blue) lines are for the 
large angle scattering cases (LAS: 0 S catt < ^ » I /71 ), and these asymp- 
totically approach the power-laws dN/dp oc p~ Te indicated by lightweight 
lines, at high and very high energies (not shown). For the 71 ft = 10 
case, also exhibited are three smaller angle scattering cases, 0 S catt < 60° 
(red), 0 S catt < 20 ° (green), and 0 S catt < 6 ° (magenta) corresponding to 
pitch angle diffusion (PAD). These have high-energy asymptotic power-law 
indices of r e = 1.65, r e = 1.99 and T e = 2.20, respectively. 

Representative particle distributions that result from our sim- 
ulation of diffusive acceleration at relativistic shocks are de- 
picted in Figure 2, highlighting several key features. These dis- 
tributions are equally applicable to electrons or ions, and so the 
mass scale is not specified. The spectral index declines and the 
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distribution is flatter for faster shocks with larger upstream flow 
(bulk) Lorentz factor 71 , when the velocity compression ratio r 
•is fixed. This is a consequence of the increased kinematic energy 
boosting occurring at relativistic shocks. Such a characteristic 
is evident, for example in the work of Kirk & Schneider (1987), 
Ballard & Heavens (1991) and Kirk et al. (2000) for the case 
of pitch angle scattering, and Ellison, Jones & Reynolds (1990), 
Baring (1999) and Ellison & Double (2004) for much larger angle 
scattering. What is much more striking in Figure 2 is that the 
slope and shape of the nonthermal particle distribution depends 
on the nature of the scattering. The asymptotic, ultrarelativis- 
tic index of T e = 2.23 is realized only in the mathematical limit 
of small (pitch) angle diffusion (PAD), where the particle mo- 
mentum is stochastically deflected on arbitrarily small angular 
(and therefore temporal) scales. In practice, PAD results when 
the maximum scattering angle 0 SC att is inferior to the Lorentz 
cone angle I/71 in the upstream region. In such cases, parti- 
cles diffuse in the region upstream of the shock only until their 
angle to the shock normal exceeds around 1 /71 . Then they 
are rapidly swept to the downstream side of the shock. The 
energy gain per shock crossing cycle is then roughly a factor 
of two, simply derived from relativistic kinematics (Gallant & 
Achterberg 1999; Baring 1999). 

To contrast these power-law cases, Figure 2 also shows our re- 
sults for large angle scattering scenarios (LAS, with 0 scat t ~ 7r ), 
where the spectrum is highly structured and much flatter on 
average than p~ 2 . The structure, which becomes extremely 
pronounced for large 71 , is kinematic in origin, where large 
angle deflections lead to the distribution of fractional energy 
gains between unity and 7 2 3 in successive shock transits by 
particles. Gains like this are kinematically analogous to pho- 
ton energy boosting by Compton scattering. Each structured 
bump or spectral segment shown in Figure 2 corresponds to an 
increment in the number of shock crossings, successively from 
1 — > 3 — > 5 — > 7 etc., as illustrated by Baring (1999); they even- 
tually smooth out to asympotically approach power-laws that 
are indicated by the lightweight lines in the Figure. The indices 
of these asymptotic results are all in the range T e < 2 . Inter- 
mediate cases are also depicted in Figure 2, with 0 sca tt ~ 4/71 . 
The spectrum is smooth, like the PAD case, but the index is 
lower than 2.23. Astrophysically, there is no reason to exclude 
such cases. From the plasma point of view, magnetic turbulence 
could easily be sufficient to effect scatterings on this intermedi- 
ate angular scale, a contention that becomes even more salient 
for ultrarelativistic shocks with 71 10 . It is also evident that 

a range of spectral indices is produced when 0 S catt is of the or- 
der of I/71 , In this case, the scattering processes corresponds 
to a transition between the PAD and LAS limits. 

Given the results of our numerical simulations, the implica- 
tions for distributions of relativistic particles in blazars are ap- 
parent. There can be a large range in the spectral indices T e of 
the particles accelerated in relativistic shocks, and these indices 
usually differ from T e ~ 2.23 . They can be much steeper, par- 
ticularly in oblique shocks (e.g., Ellison & Double 2004; Baring 
2004). However, they can also be much flatter, so that quasi- 
power-law particle distributions p~ Te with T e < 2 are read- 


ily achievable. Such flat distributions from relativistic shock 
acceleration have not usually been admitted when considering 
properties of blazar jets and their possible emission spectra. 23 

4. CONCLUSIONS 

This finding that relativistic shock acceleration produces par- 
ticle spectra with a significant range of spectral indices, includ- 
ing those with T e < 2 corresponding to inverse Compton 7-ray 
spectra with T s < 1.5, has various consequences. The consider- 
able diversity in the values T e produced in relativistic shocks 
is matched by the diversity in the intrinsic spectral indices of 
blazars indicated in Table 1. Moreover, particle distributions 
with T e < 2 are consistent with the inferred values for the three 
most distant blazars listed in the table. A hard TeV 7-ray spec- 
trum with a value T s < 2 within the context of SSC model build- 
ing (see, e.g ., Stecker, De Jager & Salamon 1996), indicates that 
the energy range of the observation is below the Compton peak 
energy in the spectral energy distribution of the source, which is 
given by E 2 times the differential photon spectrum. For extreme 
blazars, this peak can easily be at an energy above 2 TeV (de 
Jager & Stecker 2002). A simple SSC model prediction then fol- 
lows. The observation of an approximate power-law spectrum 
in the sub- TeV energy range should imply approximately the 
same index as the synchrotron emission in the optical to X-ray 
band. There are presently no high quality optical/ X-ray obser- 
vations of the distant blazars 1ES 1218+30 and 1ES 1101 -232 
simultaneous with the TeV observations. 

Our reexamination of blazar 7-ray spectra in the light of rela- 
tivistic shock acceleration theory has important implications for 
constraining the flux of the EBL in the near infrared. Specifi- 
cally, the low values of T e < 2 readily obtained in our numerical 
results implies an increase in the upper limit on the near infrared 
EBL to values above that obtained by Aharonian et al. (2006). 
Such a result is consistent with the fast galaxy evolution model 
which appears to be favored by the Spitzer observations. 

Table 1 hints at a redshift evolution of TeV blazars with a 
trend toward flatter spectra and higher isotropic luminosities at 
higher redshifts. Although only ten blazars are listed in the ta- 
ble and selection effects are important, one may speculate as to 
whether there is a general trend in blazar activity in the form of 
higher jet Doppler factors at higher redshifts. Future combined 
GLAST - TeV broadband spectral data will further define intrin- 
sic source properties, thus enabling the further investigation of 
possible redshift evolution of blazar flux and spectral character- 
istics. The H.E.S.S. and MAGIC atmospheric Cerenkov TeV 
telescopes have had remarkable success in observing blazars. 
With the VERITAS atmospheric Cherenkov telescope now on- 
line, the population of known TeV blazars will be extended con- 
siderably. 
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2 It has been suggested that electron distributions with r e < 2 can be obtained by stochastic acceleration in combination with boundary layer particle 
trapping which produces a pileup effect (Ostrowsky 2000). 

3 Our results require the implicit assumption that the electron spectra produced during blazar flares are not significantly affected by cooling by synchroton 
radiation. This requires that t acc < t cco i which constrains both the magnetic field strength and the electron energy (Baring 2002). 
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